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ABSTRACT 

Based on our recent MHD simulations, first, a concept of the successive merging of plasmoids and 
fragmentation in the current sheet in the standard flare model is presented. Then, using a 2.5-D 
electromagnetic particle-in-cell model with free boundary conditions, these processes were modelled 
on the kinetic level of plasma description. We recognized the plasmoids which mutually interacted 
and finally merged into one large plasmoid. Between interacting plasmoids further plasmoids and 
current sheets on smaller and smaller spatial scales were formed in agreement with the fragmenta- 
tion found in MHD simulations. During interactions (merging - coalescences) of the plasmoids the 
electrons were very efficiently accelerated and heated. We found that after a series of such merging 
processes the electrons in some regions reached the energies relevant for the emission in the hard 
X-ray range. Considering these energetic electrons and assuming the plasma density 10 9 -10 10 cm -3 
and the source volume as in the December 31, 2007 flare (Krucker at al. 2010), we computed the 
X-ray spectra as produced by the bremsstrahlung emission process. Comparing these spectra with 
observations, we think that these processes can explain the observed above-the- loop-top hard X-ray 
sources. Furthermore, we show that the process of a fragmentation between two merging plasmoids 
can generate the narrowband dm-spikes. Formulas for schematic fractal reconnection structures were 
derived. Finally, the results were discussed. 

Subject headings: Plasmas — Solar flares — Acceleration of particles 



1. INTRODUCTION 

It is commonly accepted that plasmoids play a very im- 
portant role in the magnetic field reconnection in solar 
flares. T heir importance, for the first time, was recog- 
nized bv lOhyama fc Shibatal ([1998ft . In the 1992 Octo- 
ber 5 flare, observed in soft X-rays by Yohkoh satellite, 
they analyzed the plasmoid which was ejected during 
the impulsive pha se upwards i n to th e corona. Study- 
ing the same flare, IKliem et al.l (|2000l ) showed that this 
plasmoid ejection was associated with the drifting pulsat- 
ing structure (DPS) on radiowaves. They proposed the 
model of this radio em ission, which was fur th er devel- 
oped in the papers by IKarlickv et al.l (l2002h. iKarlickvl 
(12004ft. iKarlickv fc Barta] (|2QQ7l ). IBarta et al.l (j2008aft . 
and IKarlickv et al.l ([20101 ). In this model, in the cur- 
rent sheet, due to tearing and coalescen ce processes 
the p l asmoi ds are formed. A s shown b y iDrake et al.l 
(|2005l. 12006ft iHoshinol (|2005ft . iPritchettl (|2006l 12008ft . 
and Karlicky (2008) during these processes electrons 
are very efficiently accelerated. The electrons are then 
trapped in plasmoids, where they generate Langmuir 
waves, which through a wave transformation produce 
the electromagnetic waves recorded on the radio spec- 
trum as DPSs. Due to limited range of plasma densi- 
ties in the plasmoid the DPS is generated in the lim- 
ited range of frequencies. In the vertical current sheet 
the plasmoids move upwards or downwards or even stay 
without any motion in d ependance on a for m of the sur- 
rounding magnetic field (Barta et al. 2008a. b). Due to a 
preference of divergent magnetic field lines in the up- 
ward direction, most of the plasmoids move upwards 
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and corresponding DPSs drift towards lower frequencies. 
Nevertheless, in some cases the plasmoids move down- 
wards and even inter act with the underlying flar e ar- 
cade as observed by iKolomanski fc Karlickvl (j2007) and 
iMillkan et all (120101) 

Recently. I Oka et al.l (|2010l ) studied the electron ac- 
celeration by multi-island coalescence processes in PIC 
model with periodic boundary conditions. They found 
that the most effective acceleration process is during the 
coalescenc e of p lasm oids (" anti-reconn ection" ) , see also 
IPritchettl (j2008h and IKarlickv fc Bartalj2007ft . 

Furthermore. IShibat a fc Tan uma (2001)) proposed that 



the current sheet, stretched by a rising magnetic rope, is 
fragmented to smaller and smaller plasmoids by the tear- 
ing mode instability in subsequently narrower and nar- 
rower current sheets (cascading reconnection) . This sug- 
gestion has been recently furt her theoretically developed 
bv lLoureiro et al.l (|2007ft and lUzdenskv et alJ(|2010f ) into 
the theory of chain plasmoid instability. An advantage 
of this concept is that it explains how very narrow cur- 
rent sheets with high current densities (requested for the 
anomalous resistivity generation and fast reconnection) 
are generated. Moreover, many X-points in this model 
give sufficient volume for an acceleration of particles. 

Besides this fragmentation describ ed by 
iShibata fc Tanumal ([20011 ). IBarta et al.l (|2010bft found a 
new fragmentation in the region between two merging 
plasmoids using MHD simulations. This fragmentation 
is caused by the tearing mode instability in the current 
sheet generated between these interacting plasmoids 
and this process is repeated in smaller and smaller 
spatial scales. This fragmentation is driven by a merging 
process of the plasmoids. 
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Considering all the above-mentioned processes, in the 
present paper, we focus our attention to two processes: 
(a) successive merging of plasmoids to large plasmoid 
and (b) fragmentation process between merging plas- 
moids. We selected these processes because we think 
that the successive merging of plasmoids can explain the 
above-the-loop-top hard X-ray source (as a large station- 
ary plasmoid). On the other hand, the fragmentation 
can explain the narrowband dm-spikes. Because both 
these phenomena are generated by accelerated electrons, 
in following simulations we use the particle- in-cell (PI C) 
model instead of MHD models (e.g. IBarta et al.ll2QQ8al lbh. 

The above-the-loop-top hard X-ray sources belong to 
the most discussed topics in recent years. The well- 
known example of such a hard X-ray sou rce is that ob- 
served in the ^30-50 keV energy range by Masuda et al. 

However, su ch events are very rare flTomczak 



12001b iPetrosian et al.ll2002t iKrucker fc : Lid l2008f ). An 
other ver y interesting e xamp le was published just re- 
cently bv IKrucker et al.l (|2010f) . They presented the hard 
X-ray source (with the energy up to ~80 keV) which 
was located 6 Mm above thermal flare loops. They de- 
rived the upper limit of the plasma density and source 
volume as n e ~ 8xl0 9 cm -3 and V ~ 8x 10 26 cm 3 , 
respectively. Just a relatively low plasma density in 
such hard X-r a y sou rces attracts attention of scientists. 
IKrucker et al.l (|201Q[ ) concluded that these hard X-ray 
sources have to be close to the acceleration region and 
the distribution function of electrons emitting hard X- 
rays is strongly non-thermal or the plasma in the source 
is very hot (up to T e ~ 200 MK). Several ideas explain- 
ing these X-ray sources have been proposed, e.g. the 
magnetic or turbulent tra pping, and dense (collisionally 
thick ) coronal sources (see iFletcherlll 9951: Ijakimiec et al.l 
119981 : IVeronig fe Brown]l2004t iPark fe Fleishman! 1201 Of ) 

The narrowband dm-spikes (further spikes) belong 
to the most interesting radio bursts due to excep- 
tionally high brightness temperatures (T5 « 10 1 5 K) 
and short durations (< 0.1 s, see the review by iBenzl 
1986). Their observational c haracte r istics were de- 
scrib e d in many papers fe.g. ISlottid Il98lt iKarlickvl 
1984; IF11 et al.l 119851: IStahli fe Magunl 119861: iBenz et al. l 
1982; IZlobec fe Karlickvl 119981 : iMeszarosova et al.l 
2003). On the other hand, the theoretical mod- 
els can be divided into two groups: a) based 
on the plasma emission and acceleration pro cesses 
(iKuiipers era!lll98iyTaiima et al.l 119901 : IWentzelll99ll : 
iBarta fc Karlicky 200 l|), and b) based on the electron- 
cyclo tr on maser (Holman~eFaL| 119801: iMelrose fc Dulkl 

IWinglee et ZT fl988: 
~M>. 



1982; Vlahos & Sharma 



1984; 



Aschwandenl 119901 : iFleishman fc Yastrebovl I1994D . To 



distinguish between these two types of models polar- 
ization and harmo nic structur e s of the spikes have 
been also studied (iGjidejl [199(1 iGiidel fc ZlobeTl fl99ll : 
IKrucker fc Benzl Il994f ). Searchi ng for a charac t eristic 
bandwidth of individual spikes Karlicky et al. (1996, 
l2QQ0h found that the Fourier transform of the dy- 
namic spectra of spikes have a power-law form with 
power- law indices close t o —5/3. Based on these results 
IBarta fc^ Karlickv (|2QQlh proposed that the spikes are 
generated in turbulent reconnect ion outflows. 

This paper is organized as follows: First, we present 
our model scenario and explain the successive merging 
and fragmentation process. Then using a 2.5-D PIC 
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Fig. 1. — Model scenario. 



model we simulate these processes. The results are then 
used in the interpretation of the above-the-loop-top hard 
X-ray sources and narrowband dm-spikes. 

2. MODEL SCENARIO AND SIMULATION MODEL 

Fig. [1] shows our model scenario, which is based on the 
'standard' CSHKP flare model (e.g. iMagara et al.|[l996l 
and references therein) . I n the central part of t he cur rent 
sheet, in agreement with Shibata fc Tanumal (|200ll ), we 
assume a fragmentation of the current sheet (stretching- 
tearing fr agmentation ). Furth ermore, based on the new 
results of IBarta et al.l (|2010bf h we propose that the re- 
connection plasma outflow (which is downward oriented) 
accumulates plasmoids in the region just above the flare 
arcade, where thus plasmoids can interact efficiently. We 
think that in some cases a large plasmoid can be formed 
here as a result of successive merging of plasmoids. On 
the other hand, between the merging plasmoids new cur- 
rent sheets are formed and in these current sheets one 
again further (but on smaller and smaller spatial scales) 
plasmoids are generated (fragmentation between merg- 
ing plasmoids). While the first process of the successive 
merging we propose for the interpretation of the above- 
the-loop-top hard X-ray sources, the second process is a 
very promising process explaining the narrowband dm- 
spikes (see the following sections I3TT1 and |3T2|) . 

For simulation we used a 2.5-D (2D3V - 2 spatial 
and 3 velocity components) fully relativistic electromag- 
netic particle-in-cell model ([Saito fc Sa kai 2004; Ka rlickvl 
l2004f ). The system size is L x x L y = 600 A x 4000 A, 
where A (=1) is a grid size. The current sheet is initi- 
ated along the line x = OA, and its half- width is L = 
10 A. In this first study we consider the neutral current 
sheet, i.e. the guiding magnetic field B z is zero. Thus, 
the initial magnetic field is: 

B = (B X1 By, B z ), 
By — —Bo for x < — L, 
By =xBo/L for — L< x< + L, 
By = Bo for x > +L, 
B x = 0,B Z =0. 
The electron-proton plasma with the proton-electron 
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Fig. 2. — The global view on magnetic field lines and corresponding current densities (reddish areas) in the x — y computational plane 
at six different times: at the initial state (a), at uo pe t = 1200 (b), at uo pe t = 2500 (c), at uo pe t = 3500 (d), at uo pe t = 4500 (e), and at uo pe t 
= 6500 (f). The x and y coordinates are expressed in A. The proton inert ial length is 40 A. 



mass ratio m p /m e = 16 is unrealistic, but taken here to 
shorten computations. Nevertheless, the electron mass 
is low enough to separate the dynamics of electrons and 
protons well. In each numerical cell located outside of 
the current sheet, we initiated no = 60 electrons and no 
= 60 protons. In the current sheet the initial number 
density was enhanced just to keep the pressure equilib- 
rium. The initial electron temperature was taken to be 
the same in the whole numerical box as T = 10 MK 
and the temperature of protons was chosen the same as 
electrons. The plasma frequency is u pe At = 0.05 (At 
is the time step which equals to 1), the electron Debye 
length is = 0.6 A, and the electron and proton in- 
ert ial lengths are d e = 10 A and di = 40 A, respec- 
tively. To study successive coalescence processes among 
several plasmoids, we initiated a formation of 10 plas- 
moids along the current sheet by a cosine perturbation 
of the electric current density in the sheet; with the k- 
vector k = 2tt- 10/4000 = 0.0157 A" 1 and the amplitude 
corresponding to the current density j given by the mag- 
netic field in the current sheet (j = V X B). We made 
computations with several initial values of the plasma /? 
parameter. Here, we present the results with f3 = 0.07. It 
gives the ratio of the plasma densities in the centrum and 
out of the current sheet as 15.3. The free boundary condi- 
tions were used. All computations were performed on the 
parallel computer OCAS (Ondrejov Clust er for Astro - 
physical Simulations), see |http: / / wave.asu.cas.cz / ocas| 

3. RESULTS 

A global evolution of the magnetic field lines and the 
corresponding electric current densities (reddish areas) in 
the system is shown in Fig. [2j When 10 small plasmoids 
were formed, at about uo pe t = 1200 (Fig. [2}d), then the 
plasmoids started to interact and merge into larger plas- 



moids. Due to free boundaries used in the system and 
small asymmetries in the initiation, the plasmoids suc- 
cessively merged into one large plasmoid formed in the 
bottom part of the system (Fig.(2f). On the other hand, 
Fig. [3] presents a more detailed view on this evolution as 
well as the distribution of numerical electrons (points), 
having the energy greater than 40 keV, at four different 
times: at uo pe t = 3200 (a), at uo pe t = 3500 (b), at u pe t = 
5600 (c), and at u pe t = 7800 (d). As can be seen here, 
in each coalescence process electrons are very efficiently 
accelerated - see an increase of the number of numerical 
electrons in dependence on time. In Fig. (3Jd and [3]i, we 
selected the region (see white crosses = the centra of the 
circles of the radius 70 A) with an enhanced number of 
accelerated electrons. In these regions we determined the 
normalized (the maximum equals to 1) electron distribu- 
tion functions in all three coordinates (Fig. 2^,b). For 
comparison in Fig. H^i the thermal distribution function 
with the temperature 66 MK is added. It shows that 
this distribution function has clear nonthermal tails. On 
the other hand, the computed distribution function in 
Fig. HJd is nearly the thermal one, see the fit of this func- 
tion by the thick full line corresponding to the thermal 
distribution function with the temperature 107 MK. 

In summary, Figs. [3] and 2] show that the successive 
merging (coalescences) of the plasmoids step by step 
increase the energy (and the temperature) of acceler- 
ated electrons. We found that at some regions and for 
short times during the coalescence process the distribu- 
tion functions deviate from the thermal ones (the power- 
law tails in Fig. U K or even bumps - on-ta il, see e.g. Fig. 4 
in the paper bv lKarlickv fc Bartal ([20071 )), but very soon 
these distribution functions are thermalized (i.e. changed 
to the Maxwellian ones) by fast wave-particle processes 
(anomalous collisions). 
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Fig. 3. — The detailed view on magnetic field lines in the x — y computational plane and distribution of numerical electrons (points) 
having the energy greater than 40 keV at four different times: at uo pe t = 3200 (a), at uo pe t = 3500 (b), at uo pe t = 5600 (c), and at oo pe t = 
7800 (d). The white crosses, located at b) x = 20, y = 900, and at d) x = 20, y = 150, show the regions where the distribution functions 
and the X-ray spectra were computed (see Figs. [4] and [5]). The x and y coordinates are expressed in A. 



3.1. Successive merging of plasmoids and the 
above-loop-top hard X-ray sources 

We computed the hard X-ray spectra generated at 
two selected regions (the white crosses in Fig. [3]) and 
compared them wi th that observ ed during the Decem- 
ber 31, 2007 flare (IKrucker et al. 1201 Oft . In accordance 
with the paper by IKrucker et al. (|2Q10h we have taken 
the plasma density in the range 10 9 — 10 10 cm -3 and 
the source volume V = 8x 10 26 cm 3 . For the com- 
putation of the hard X-ray spectra we used two meth- 
ods: (A) the non-thermal bremsstrahlung method, in 



which the spectrum was computed as a sum of contri- 
butions of the bremsstrahlung emission of all numeri- 
cal electrons, fo r details see the rel a tions (10) and (11) 
in the paper by iKarlickv fc Kosugil (|2004f ). and (B) the 
thermal b remsstrahlung method for specific p lasma tem- 
peratures (|Tandberg-Hanssen fc Emslielll988t ) . The com- 
puted spectra together with the observed spectrum are 
in Fig. [5ji, b. The spectra computed by the method (A) 
and (B) (in Fig. (5Jd) are similar. The observed spec- 
trum is between the spectrum computed for the source 
plasma density 10 9 cm -3 and that for 10 10 cm -3 , espe- 
cially in later phases of the model evolution. But, the 
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Fig. 4. — The normalized electron distribution functions (thin 
full line means f(v x ), dashed line f(v y ), and dotted line f(v 2 )) com- 
puted in the regions shown by the white crosses in Fig. [3] at two 
different times: at uj pe t = 3500 (a) and at (jj pe t = 7800 (b). The 
thick full line in (a) means the thermal distribution function with 
the temperature 66 MK and the thick full line in (b) the distribu- 
tion function with the temperature 107 MK, which roughly fits the 
computed distribution functions. 



Fig. 5. — The X-ray spectra (a,b) (the full line for the source den- 
sity n e = 10 9 cm -3 , and the dashed line for n e = 10 10 cm -3 ) cor- 
responding to the distributions in Fig.[J^i, b, computed by method 
(A). The X-ray spectra in (b) (the dashed-dot line for the source 
density n e = 10 9 cm -3 , and the dashed-dot-dot-dot line for n e = 
10 10 cm -3 ) corresponding to the distribution in Fig.^J), but com- 
puted by the method (B) for the temperature 107 MK. For com- 
parison the X-ray spectrum (dotted line) observed during the De- 
cember 31, 2007 flare (according to Krucker et al. 2010) is added. 



observed spectrum has a different form, comparing to 
the computed ones. Considering our results we think the 
observed X-ray power-law spectrum is given by a sum 
of emissions from many locations with different thermal 
and nonthermal distribution functions. 

3.2. Fragmentation between merging plasmoids and 
narrowband dm-spikes 

In accor d ance w ith the results presented by 
Bart a et al.l jjjlOb), we studied the structure of 
the magnetic field in the region between merging large 
plasmoids. One example of such a structure, formed 
at the time u pe t — 7200 is presented in Fig. [6l The 
reddish areas mean the electric current densities with 
the initial orientation of the electric current, and the 
blue ones mean those with the opposite orientation; red 
and green lines are positions of B y = and B x = 0, 
respectively; their intersections represent the X- and 
O- type null points. The figure shows that the initial 
current sheet is fragmented into several sub-current 
sheets. In these secondary current sheets further tearing 
can exist provided they are sufficiently long. Our results 
thus indicate, that the fragmentat ion cascade seen 
in a large-scale MHD simulation (|Barta et al.l l2010bl ) 
continues actually down to the dissipation scale of the 
order « d{ (the thickness of the CS at the bottom panel 
of Fig. [6] is ~ 30 A). At this scale, however no further 
fragmentation has been observed, the current structures 
dissipate by plasma kinetic processes 'silently'. 

Considering the structures of fragmented current 
sheets (cascade of plasmoids) obtained by numerical sim- 
ulations (Fig. [6]), we constructed the schematic struc- 
ture of fragmentation. An example of such a structure is 
shown in Fig. [71 The circles mean plasmoids with posi- 
tively (+) and negatively (— ) oriented electric currents, 
the parallel lines mean boundaries of current sheets. In 
such a structure the radius of the plasmoid Ri can be 
written as: 

Ri = AR i+ i + BR i+ 2, (1) 

where A and B are the numbers of plasmoids and current 
sheets in specific current sheet, i means the index of the 
plasmoids with the same size. The number of plasmoids 



increases with i as ~ B\ Then, the ratio between 
subsequent radii of plasmoids in the plasmoid cascade 
can be expressed by the infinite continued fraction: 



Ri 
Ri+l 



B 



A 



(2) 



Now, assuming that the energy in the plasmoids is pro- 
portional to their area, the dependence of the function 
Ei = Rf x Hi on the k- vector scale (ki = 27r/Ri) can be 
computed. This function is the power-law one with the 
power- law index 



P : 



log(C • B) 
log 



2, 



(3) 



where C is the free parameter expressing possible devi- 
ations from our assumption about the plasmoid energy. 
For the structure presented in Fig. [71 and for C = 1 this 
power-law index is p = —1.54. 

The power-law dependence were found also in the anal- 
ysis of t he frequency bandwidth o f the narrowband dm- 
spikes ([Karlickv et al.lfl996L [2000). Therefore, consider- 
ing this fact, the model scenario show n in Fig. [H and the 
turbu lent-plasma model of spikes by iBarta fc Karlickvl 
(|200lh . we propose that the narrowband dm-spikes are 
generated in fragmentation processes between merging 
plasmoids in the reconnection outfl ow in the region 
above the flare loop arcade. While Bar ta fc Karlickvl 
(|200lf ) supposed (silently) MHD/ Alfvenic turbulence in 
the super-sonic outflows, our recent simulations indi- 
cate that the tearing/coalescence cascade might be more 
likely source of turbulence. However, in a fully developed 
MHD turbulence all plasma wave- modes are present, 
anyway. In the fragmentation region, the plasmoids of 
all spatial scales can interact and accelerate electrons. 
These electrons are trapped in plasmoids of different 
sizes. In each plasmoid they can generate the radio 
emission in the frequency range corresponding to the 
range of the plasma density in this plasmoid. Due to 
expected power-law dependence of spatial scales of these 
fragmented plasmoids, the dependence of bandwidths of 
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Fig. 6. — Fragmentation of the current sheet in the direction 
perpendicular to the initial current layer, i.e. in the current sheet 
between two large interacting plasmoids, at uj pe t = 7200. The 
reddish areas mean the electric current densities with the initial 
orientation of the electric current, and the blue ones mean those 
with the opposite orientation; red and green lines are positions of 
By = and B x = 0, respectively; their intersections represent the 
X- and O- type null points. 



the resulting radio emission should be the power-law one 
as observed in the narrowband dm-spikes. 

To support this idea, in Fig. we present the radio 
spectrum observed during the 28 March 2001 flare by 
two O ndfejov radiospectr ographs (0.8 — 2.0 and 2.0 — 4.5 
GHz) (Jifick aet al.l fl993). It shows two positively drift- 
ing pulsating structures (DPSs) which according to our 
previous studies indicate two plasmoids moving down- 
wards to the sources of the narrowband dm-spikes (gen- 
erated in the region of fragmentation, compare to the 
model scenario in Fig. [1]). The mean Fourier spectrum 
of these spikes is the power-law one with the power-law 
index —1.5. 

4. DISCUSSION AND CONCLUSIONS 



x 



Fig. 7. — Example of the fractal reconnect ion structure corre- 
sponding to the relations JT} and (0, for A=2 and B=3. Double 
straight lines delineate current sheets with interacting plasmoids 
(circles). The plus and minus signs express orientation of electric 
currents in the plasmoids. X and Y are in arbitrary units. 



Important aspect of our model is that we used free 
boundary conditions which enabled successive merging 
of small plasmoids into a large final plasmoid. In sim- 
ulations, we recognized also a fragmentation of current 
sheet between two merging plasmoids. 

We showed that these processes very efficiently accel- 
erate electrons to energies relevant for the emission in 
the hard X-ray range. Based on this result we propose 
that the above-the-loop-top hard X-ray sources are pro- 
duced by the successive merging of plasmoids in the re- 
gion with the turbulent reconnection outflow, in the re- 
gion just above the loop arcade. Computed X-ray spec- 
tra support this idea. To explain a difference between the 
slopes of the computed and observed X-ray spectra we 
propose that the observed power-law spectrum is a sum 
of emissions from many locations with different thermal 
and nonthermal distribution functions. 

In simulations we used the particle-in-cell model. Al- 
though, the studied processes are self-similar (i.e. they 
does not depend on scales, see e.g. the MHD and PIC 
simulations in the paper bv lTajima et al.l ([1987)), the re- 
sults need to be taken with a caution, because the spatial 
and time scales in the model and real plasmoids differ in 
several orders of magnitude. On the other hand, it is 
beyond a possibility of any present numerical model to 
take into account all these scales. 

Further aspect of the PIC modelling is that the 
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Fig. 8. — The radio spectrum observed during the 28 March 
2001 flare by two Ondfejov radiospectrographs (0.8-2.0 and 2.0- 
4.5 GHz) supporting our model. It shows the drifting pulsat- 
ing structures (DPSs) which drift towards narrowband dm-spikes. 
The 0.8-2.0 GHz spectrum was shortly interrupted at 12:08:34- 
12:08:36 UT. 



Coulomb collisions were not considered in our model. 
Namely, the time interval of our computations is much 
shorter than the collision time in the-above-the-loop-top 
X-ray sources; e.g. u pe t = 7800 for the plasma den- 
sity n e = 10 9 cm -3 corresponds to 4.5 x 10 -6 s and 

the collision time to 0.115 T 3 7 /2 s, where T 7 = T e /10 7 . 
On the other hand, the Coulomb collisions are essen- 
tial for the bremsstrahlung X-ray emission. The hard 
X-ray spectra in the December 31, 2007 flare were de- 
tected as the mean ones over intervals of several seconds. 
Furthermore, the-above-the-loop-top hard X-ray source 
lasted of about 2 minutes. Thus comparing these times 
with the collision one, the collisions influence not only 
the observed spectra, but also a duration of the X-ray 
source. Therefore some re- acceleration of electrons is 
even needed. In the present model such a re-acceleration 
is very probable, because the plasmoids are generated in 
a broad range of spatial scales and electrons can travel 
several times though regions of interacting plasmoids. 
Although an inclusion of the Coulomb collisions and cor- 
responding prolongation of computations are beyond the 
possibility of our present PIC modelling, we think that 
the observed spectra are given by a competition of the 
collisions with the re-acceleration of electrons. On the 
other hand, in PIC model the anomalous collisions (wave- 
particle interactions), which are much more effective than 
the Coulomb collisions, are present as confirmed by fast 
thermalization of nonthermal distribution functions. 

We made additional computations with different initial 
parameters and we found: a) the energy gain of acceler- 
ated electrons increases with the decrease of the plasma 
beta parameter, and b) the increase of the proton- 
electron mass ratio m p /m e makes computations longer, 
but results are similar. 

We compared the present simulation also with that in 
the numerical model which size was two times smaller 



(L x x L y = 600A x 2000A) and in which only 5 plas- 
moids were initiated (contrary to 10 plasmoids in the 
present simulation). In this case the final mean energy 
of accelerated electrons was 5.3 times greater than the 
initial one, compare with that of 10.7 times (from the 
initial temperature 10 MK to final temperature 107 MK) 
in the present case. Namely, each coalescence process in- 
creases the energy of accelerated electrons, therefore the 
number of successive coalescence processes is essential for 
their final energy. 

For calculations of the hard X-ray spectra, presented 
in Fig. [5b we used two methods. The obtained results 
are similar. Small differences are due to differences in 
these methods and deviations of the computed distribu- 
tion functions from the thermal one. 

The plasmoids in 2-D are in reality 3-D magnetic ropes. 
While in 2-D the trapping of energetic electrons is a 
natural consequence of a close magnetic field structure 
of the plasmoid, in 3-D, this structure is only semi- 
closed. However, we consider the merging processes in 
the turbulent reconnection outflow therefore the mag- 
netic trapping of elec trons, similar to that proposed by 
iJakimiec et al I Jl998) is highly probable. Moreover, the 
coalescence fragmentation process, which generates the 
reverse electric currents (which in 3-D has to be closed 
in finite volume) will contribute to a full trapping of elec- 
trons. 

In a greement with the conclusions by iKrucker et al.l 
(2010), in the model the acceleration region is very close 
to the hard X-ray source. It enables to re-accelerate en- 
ergetic electrons, which loss their energy due to colli- 
sions. Acceleration regions are among interacting plas- 
moids and also between the plasmoids and the arcade of 
flaring loops. This model can explain not only the above- 
the-loop-top hard X-ray sources, but also the loop-top 
sources because the arcade of loops is, in principle, the 
'plasmoid' fixed in its half height at the photosphere. 

Considering all aspects of the fragmentation process 
(power-law spatial scales of plasmoids, effective acceler- 
ation of electrons, trapping of electrons in plasmoids, lo- 
cation in the reconnection plasma outflow) we think that 
this process can explain a generation of the narrowband 
dm-spikes. We supported this idea by the radio spectrum 
observed during the 28 March 2001 showing DPSs drift- 
ing towards the narrowband dm-spikes. Furthermore, 
it is known that more than 70 % of all groups of dm- 
spikes are observed during the GOES-rising-flare phases 
([Jificka et al.ll200T| ). Although these arguments support 
the presented idea, further analysis of the narrowband 
dm-spikes and their modelling is necessary. 

In this first study, we considered only the neutral cur- 
rent sheet, i.e. B z = 0. For more realistic description, 
we plan to extend our study also to cases with non-zero 
guiding magnetic field. 

The presented model is a natural extension of our pre- 
vious models explaining the plasmoid formation, its ejec- 
tion, and corresponding DPS. The question arises why 
the above-the-loop-top hard X-ray sources are very rare 
comparing, e.g., with DPSs or dm-spikes observations. 
We think that the above-the-loop-top hard X-ray source 
is large and stationary plasmoid, which is sufficiently 
dense and in which there is sufficient amount of ener- 
getic electrons. Positional stationarity and location of 
this plasmoid is given by surrounding magnetic field and 
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location whe re this plasmoid sta rts to be formed, see 
the paper bv iBarta et al.l (|200 8b). On the other hand, 
the plasmoids generating DPSs or dm-spikes need not to 
be so dense and they need not to have such amount of 
energetic electrons. It is known that for generation of 
the radio emission a number of energetic electrons can 
be much smaller than that for the hard X-ray emission. 
Namely, the intensity of the radio emission depends on 
derivatives of the electron distribution function in the 



momentum space; not on the absolute amount of ener- 
getic electrons as in the case of the hard X-ray emission. 
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of Astronomical Institute of the Czech Academy of Sci- 
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